AD-A072  199 


UNCLASSIFIED 


ARMY  MATERIALS  AND  MECHANICS  RESEARCH  CENTER  WATERTOWN  MA  F/G  20/11 
ILLUSTRATIONS  OF  TWO  SIMPLE  METHODS  OF  DETERMINING  STRESS  INTEN— ETC(U) 
APR  79  F I BARATTA 

AMMRC -MS-79-3  NL 


1 


DDC  _F1LE  JflPY.  SB  A 0721  99 


AMMRC  MS  79-3 


nanEfn 


M 


Approved  for  public  release;  distribution  unlimited. 


ILLUSTRATIONS  OF  TWO  SIMPLE  METHODS  OF 
DETERMINING  STRESS  INTENSITY  FACTORS  FROM 

(a)  STRESS  CONCENTRATIONS  AND 

........  - ' ■ " * 

lb)  STRESS  DISTRIBUTIONS 


April  1979 


ARMY  MATERIALS  AND  MECHANICS  RESEARCH  CENTER 
Watertown,  Massachusetts  02172 


79  08  1 022 


MNkAMSI 


f 


A 

* 


The  findings  in  this  report  are  not  to  be  construed  as  an  official 
Department  of  the  Army  position,  unless  so  designated  by  other 
authorized  documents. 

Mention  of  any  trade  names  or  manufacturers  in  this  report 
shall  not  be  construed  as  advertising  nor  as  an  official 
indorsement  or  approval  of  such  products  or  companies  by 
the  United  States  Government. 


DISPOSITION  INSTRUCTIONS 


* 


Destroy  this  report  when  it  it  no  longer  needed. 
Do  not  return  it  to  the  originetor. 


MTV  ClAMlPlCATtON  OP  TMII  P*01(Ml«n  9a M 


Block  No.  20 


ABSTRACT 


Illustrative  examples  of  two  methods  of  determining  approximate 
stress  intensity  factors  from  (a)  stress  concentration  factors  and 
(b)  from  stress  distributions  are  presented. 

In  the  first  case  several  example  problems  are  presented  whereby 
stress  intensity  factors  are  obtained  from  known  stress  concentration 
factor  data.  These  illustrative  problems  include  data  obtained  from 
both  closed  form  and  experimental  solutions.  The  results  of  this  method 
are  compared  to  a well-known,  accurate  analysis,  and  guidelines  regard- 
ing such  usage  are  presented. 

In  the  second  case  several  approximate  methods  of  determining 
stress  intensity  factors  utilizing  stress  distributions  available  in  the 
literature  are  presented  to  the  reader.  Example  problems  are  illustrated 
and  compared  to  more  accurate  results  when  available.  Estimates  of  the 
accuracy  and  limited  guidance  on  the  usage  of  such  methods  are  given. 
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FOREWORD 


The  two  illustrative  methods  of  determining  stress  intensity  factors  de- 
scribed herein  were  previously  written  as  part  of  a chapter  to  be  included  in 
an  AGARDagraph  report  entitled  "On  Practical  Applications  of  Fracture  Mechanics." 
Since  these  two  methods  are  of  interest  to  the  technical  community  it  was  con- 
sidered worthwhile  to  publish  them  as  a monograph. 


It  is  also  intended  that  the  results  summarized  in  the  present  document 
will  become  part  of  a more  complete  standardization  type  handbook  that  will 
provide  easy  access  to  stress  intensity  factors  directly  related  to  Army 
applications. 
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PART  A.  STRESS  INTENSITY  FACTORS  FROM  STRESS  CONCENTRATIONS 


INTRODUCTION 

Stress  intensity  factors  K may  be  computed  from  stress  concentration  factors 
as  the  notch  radius  p approaches  zero.  This  may  be  observed  by  noting  that  a crack 
with  a small  but  finite  tip  radius  p still  has  its  tip  embedded  within  the  stress 
field  as  described  by  the  elastic  stress  field  equations  for  mode  I given  below1 2 
as  the  crack  tip  is  approached,  i.e.,  as  r 0: 

Ki  6,,  .0  . 30. 

oxx  = ■ A cosj(l-siny  siny) , 

/2itr 

Ki  0 q 70 

Oyy  = — cosyO-esin-^-  sin^— ),  and  (1) 

/2irr 

Kt  . 0 0 30 

T_v  = sin?  COST-  COST—. 

xy  2 2 2 

Of  course,  very  near  the  tip  itself,  the  finite  radius  will  alter  the  stress  dis- 
tribution. However,  the  maximum  stress  omaX  will  occur  on  the  notch  and  must  be 
proportional  in  intensity  to  the  stress  field  around  it,  i.e.: 

°max  * KI-  (2) 

The  only  characteristic  dimension  available  to  make  Equation  2 complete  is 
the  notch  radius  p,  consequently 

KI  * Q°max  ^ (3) 

where  Q is  a numerical  constant. 

For  an  elliptical  hole  stressed  by  ag  in  an  infinite  sheet,  the  stress  concen- 
tration equation*  is: 


-Jr  - [l+2(c/ p)^] , (4] 

where  2c  is  the  major  axis  of  the  ellipse. 

The  stress  intensity  equation  for  a through  crack  of  length  2a  which  the  el- 
lipse degenerates  to  in  the  limit  is: 

Kj  = o0  Cna)*.  (5) 


1.  PARIS,  P.  C..  »nd  SIH,  C.  C.  Stress  Analysis  of  Cracks.  A STM  Special  Technical  Publication  381,  1970.  p.  30-80 

2.  FNGUS,  C.  E.  Stresses  in  a /late  Hue  to  the  Presence  of  Cracks  and  Sharp  Comers.  Tran*.  Instn.  Nav.  Archil.,  v.  55.  1913,  p.  219. 


Substitution  of  Equations  4 and  S into  Equation  3 and  taking  the  limit  as 
p -*•  0,  gives  Q ■ j— , as  indicated  in  the  following: 

S H 

°0 (*•)**  ■ Q oo[l+2(c/p)**]  p*1  or  Q ■ VC-^  v “ 5“ 

p’+2c’ 

and  therefore  Equation  3 becomes,  as  deduced  by  Erwin:3 


Kj  - lim  kt  pop*1. 
p 0 I 


for  mode  I only,  provided  Kjj  ■ Kjjj  ■ 0;  and  kt  is  the  stress  concentration 
factor,  i.e.,  k,  ■ o /oo. 

I ‘I 

Formulas  corresponding  to  Equation  6 can  be  derived1 for  Kjj  for  Km  as 
well  as  Kg,  which  is  a special  case  of  mode  I,  and  these  are: 

KtI  ■ lim  r/H  kt  t op^ , (' 

p -►  0 II 

provided  K.  ■ KTII  - 0,  where  kt  is  the  in-plane  shear  stress  concentration 

II 

factor,  i.e.,  kt  - x-~* . 

rII  T0 


Km  ■ lim  /it  kt  t0pV 
p 0 III 


provided  Kj  * Kjj  ■ 0,  where  kt^  3S  °ut*°f-pIBn*  shear  stress  concentration 


factor,  i.s.,  kt 


III  T0 


Finally, 


Kfi  - l/2((3*v)/(l*v))  lim  Hwix 

p •+■  0 

where  v is  Poisson's  ratio  and  "max  is  the  maximum  bending  moment. 


3.  IRWIN,  G.  R.  Future  Mechanics.  Structural  Mechanics.  Goodtfr  ami  Hoff,  ml.,  Prrgamon  Ptpm,  Oxford,  IdftO, 


Equations  6 to  9 provide  powerful  but  simple  tools  to  obtain  stress  intensity 
factors  directly  from  stress  concentration  factors,  because  included  in  these  limit 
equations  are  stress  concentration  factors.  Thus,  if  the  stress  concentration  fac- 
tor as  a function  of  p for  a given  loading  system  is  known,  by  taking  p to  the  limit 
in  the  appropriate  equation  given  above,  the  stress  intensity  factor  can  be  obtained 

There  are  a number  of  references  available,  such  as  Heywood,1*  Neuber,4 5 
Peterson6  and  others,  that  provide  experimental  stress  concentration  factor  data  as 
well  as  theoretical  expressions  for  many  configurations.  The  data  taken  from  some 
of  these  references  as  well  as  others  will  be  utilized  in  the  examples  to  follow. 


EXAMPLE  PROBLEMS  AND  APPLICATIONS 

To  illustrate  the  approach,  consider  the  semi-elliptical  edge  notch  of  depth  c 
as  stressed  by  og  in  a semi-infinite  sheet  shown  in  Figure  1.  The  stress  concen- 
tration factor  k from  Bowie7  is  given  in  Table  1.  By  reworking  the  data  from  Refer- 
ence 7 (third  and  fourth  columns)  so  that  it  is  dimensionless,  and  plotting  kt/p/c 
as  a function  of  p/c  as  shown  in  Figure  2 (a  semi -log  plot),  we  can  graphically 
determine  the  limit  on  Kj  as  p 0.  It  is  seen  from  the  figure  that  as  p ■*•(), 
c -*•  a,  where  c is  the  crack  length.  The  extrapolated  value  of  ao/2) 

approaches  2.26  and  thus. 


*1 

/tuT  Oq 


2.26 

2 


1.13. 


Figure  1 . Semi-elliptical  edge  notch  in  a semi-infinite  sheet. 


4.  HEYWOOD,  R.  B.  Designing  by  Photoelasticity.  London,  Chapman,  and  Hall,  Ltd.,  1952. 

5.  NEUBER,  H.  Theory  of  Notch  Stress:  Principles  for  Exact  Stress  Calculation  of  Strength  with  Reference  to  Structural  Form 
and  Materials.  AEC  Technical  Report,  TR  No.  4547,  1958. 

6.  PETERSON,  R.  E.  Stress  Concentration  Factors.  John  V/iley  A Sons,  New  York,  1974. 

7.  BOWIE',  O.  L.  Analysis  of  Edge  Notches  in  a Semi-Infinite  Region.  S.  Math,  and  Physics,  v.  45,  no.  4,  December  1966;  also 
Army  Materials  and  Mechanics  Research  Center,  AMRA  TR  66-7,  June  1966. 


Table  1.  REDUCTION  OF  STRESS  CONCENTRATION 
FACTOR  TO  STRESS  INTENSITY  FACTOR 


Figure  2.  Stress  intensity 
factor  from  theoretically 
determined  stress  concert 
tration  factors. 


The  more  accurate  value8 9 10  is  1.1J15  and  the  error  is  less  than  •*■1%.  If  the  stress 
concentration  is  known  for  extremely  sharp  notches  as  accomplished  in  Reference  9 
extrapolation  is  eliminated  and  the  error  can  be  less  than  1\ . 


Experimental  data  can  be  handled  in  a similar  fashion  as  that  given  in  the 
preceding  example.  Consider  the  experimental  data  given  by  Frocht1^  for  relatively 
deep  symmetrical  U-notches  in  a finite  plate  of  notch  depth-to-plate  width  L/b  of 
1/2.  This  data  is  now  presented  in  dimensionless  form  in  Figure  3;  /p/L  k is 
plotted  as  a function  of  p/L  on  semi -log  paper. 


8.  HOTTER.  W.  T.  Oh  the  flexural  RigiJitx  of  a Ream  h'eakenej  by  7>mjiw  Saw-thtts.  I.  Appl.  Mfoh..  v.  32.  I'M.  5,  p,  237, 

9,  HASEBE.  N„  and  KITANDA,  Y.  Calculation  of  Stress  Intensity  Fecit*  pom  Stress  Ctmemtralum  Factor.  Eng.  Fracture  Meeh., 
v.  10.  1978.  p.  215-221. 

10.  FROCFTT.  M.  M.  4 Photoeiestie  Inrestifahoa  of  Stress  Cimcentnatiims  Hue  to  Small  Fillets  ami  Or\u»ys  in  Tension.  Tech.  Notes 
Nat.  Advti.  Comm.  Aeronaut.,  Washington,  Report  No  2442,  1951. 
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k&t 

T.O 

3.065 

1.0 
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4.63 

.333 
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9.0 

7.41 

.111 
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16.0 

9.625 

.063 

2.41 

32.1 

13.320 

.031 

2.35 

SI  .02 

20.77 

.012 

2.28 

361.4 

43.2 

.0028 

2.29 

Figure  3.  Stress  intensity 
factor  from  experimental 
stress  concentration  data. 


Obviously  the  data  point  at  p/l.  = 0.030  is  an  "outlier";  particularly  since 
the  notch  acuity  is  quite  sharp,  which  can  result  in  a large  experimental  error 
when  attempting  to  obtain  the  stress  concentration  factor  by  photoelasticity. 
Therefore,  this  data  point  was  ignored  in  Figure  3 and  a straight  line  was  drawn 
through  the  remaining  data  resulting  in  k/oo*'rTa  = 1.19  (note  as  p -*■  0 , L -*■  a") . 

The  more  accurate  value  is  1.15  from  Reference  11  resulting  in  an  error  of  + 3V 

If  all  of  the  data  in  Figure  3 had  been  utilized  to  determine  kj  by  the  above 
extrapolation  technique,  considerable  error  would  have  resulted.  By  judiciously 
eliminating  the  data  point  (s')  that  appear  to  be  in  error  a reasonable  approxima- 
tion to  Kj  can  be  realized. 

Finally  to  illustrate  how  Equation  o can  be  applied  to  a closed-form  expres- 
sion, consider  the  formula  given  in  Reference  12  which  describes  the  stress  con- 
centration factor  for  symmetrical  side-edge  U-notches  stressed  at  infinity  by 
oq , as  shown  in  Figure  3: 


k,  - [0. 780+2. 243/U7p] [0 .993*0. 180(L/b) -1 .060(L/b) 2+l .710(L/b) 3]  (10) 

where  1.0  $ L/p  < 361  and  0 < L/b  S 0.723. 


II.  BOWIE,  O.  L.  RfftanpiUr  Ttnsilt  Shftt  with  Symmetric  K/jrr  (Vr'Ij.  Trans.  ASME, J.  Applied  Meek.,  v.  31,  set.  E,  no.  2, 

June  1969. 

I }.  BARATTA,  F.  I..  and  NEAL,  P.  N.  .Vfrm-CtwtcenfmMrwi  h'tchvt  in  I Shj/H-,1  jnJ  SrmiMlipt1c*l  Stuped  Kfjre  S'otckn.  J.  Strain 
Analysis,  v.  S.  no.  2.  1970,  p.  121-127;  also  Army  Materials  ami  Mechanics  Research  ('enter.  AMMRO  TR  *1-1,  January  1970. 
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A comparison  of  the  results  given  by  liquation  12  to  the  more  exact  data 
given  in  Reference  11  is  shown  in  Table  2.  It  is  seen  that  the  error  in  compar- 
ing Kj/'/rra’  o0  from  Equation  12  to  that  of  Reference  11  is  no  greater  than  +1%. 

It  is  concluded  that  the  method  of  obtaining  stress  intensity  factors  from 
stress  concentration  data  is  simple  to  employ  but  the  accuracy  is  dependent  upon 
the  accurate  determination  of  the  stress  concentration  factor  as  a function  of 
the  notch  or  concentrator  acuity.  Therefore,  the  user  must  provide  judgment  to 
obtain  reasonable  results  when  utilizing  this  technique. 


Table  2.  COMPARISON  OF  STRtsS  INTENSITY 
FACTOR  U[/*-a  j0)  DATA 


a/b 

From 
tq.  12 

From 
Ref.  11 

Differ- 
ence (t) 

0.112 

1.12 

1.12 

0 

.1667 

1.12 

1.12 

0 

.2222 

1.12 

1.13 

0 

.3333 

112 

1.13 

0 

.3889 

1.13 

1.14 

+ 1 

.5001 

1.16 

1.15 

+1 

.7230 

1.36 

1.37 

+ ] 

SUMMARY 


The  method  of  obtaining  Kr  from  kr  data  is  simple  and  powerful.  But  accu- 

1 TI 

racy  is  dependent  upon  the  accurate  determination,  either  experimentally  or  ana- 
lytically obtained,  of  the  stress  concentration  factor  as  a function  of  notch 
acuity.  The  accuracy  of  those  cases  examined  was  well  within  that  required  for 
engineering  applications.  However,  one  should  know  "a  priori"  how  kt  generally 

behaves  as  a function  of  notch  radius  to  notch  depth  and  eliminate  or  include 
appropriate  values  for  a realistic  determination  of  stress  intensity  factor  at 
the  limit  of  p/c. 
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PART  B.  STRESS  INTENSITY  FACTORS  FROM  STRESS  DISTRIBUTIONS 


INTRODUCTION 

There  are  several  methods  by  which  stress  intensity  factors  are  estimated  from 
stress  distributions.  For  example,  Williams  and  Isherwood1  approximate  the  strain 
energy  release  rate  in  terms  of  a mean  stress  attributed  to  the  appropriate  geometry 
discontinuity  for  various  problems  of  interest.  Kobayashi2  applies  the  stress  dis- 
tribution at  the  crack  face  through  the  use  of  Green's  function  in  conjunction  with 
the  elimination  of  residual  surface  tractions  which  arise  because  of  the  geometric 
discontinuity  on  the  crack  surfaces  located  in  the  uncracked  region.  Shah3  applies 
a similar  procedure  but  in  addition  uses  the  superposition  method.  Cartwright,4 5 
rather  than  applying  the  stress  distribution  to  the  crack  face,  has  estimated 
stress  intensity  factors  by  applying  the  stress  due  to  the  geometry  discontinuity 
to  only  the  crack  tip. 

These  approaches  will  be  demonstrated  in  the  following  section  by  solving 
several  problems  and  comparing  these  results  to  more  exact  analyses  where  available. 
(Cracks  inclined  to  the  applied  stresses  are  not  considered  here.) 

METHODS  AND  APPLICATIONS 

Williams  and  Isherwood1  determine  the  strain  energy  release  rate  from 

G = ircrm2a/E  (1) 

where  am  is  the  mean  stress,  a is  the  crack  length,"  and  E is  Young's  modulus,  by 
application  of  the  mean  stress  am  over  the  crack  length,  assuming  no  crack  is 
present.  This  approach  can  be  illustrated  by  considering  Figure  la,  an  infinite 
plate  stressed  by  uniaxial  tension  oq*  containing  two  radial  cracks  at  the  bound- 
ary of  the  internal  circular  hole.  The  mean  stress  for  this  geometric  disconti- 
nuity is 


where  from  Reference  5 the  stress  distribution  a(r)  is: 

a(r)  = o0 [l+(l/2) (R/r)2+(3/2) (R/r)4] . (3) 


1.  WILLIAMS,  J.  G.,  and  ISHERWOOD,  D.  D.  Calculation  of  the  Strain-Energy  Release  Rates  of  Cracked  Plates  by  an  Approximate 
Method.  J.  Strain  Anal.,  v.  3,  no.  1,  1968,  p.  17-22. 

2.  KOBAYASHI,  A.  S.  A Simple  Procedure  for  Estimating  Stress  Intensity  Factor  in  Region  of  High  Stress  Gradient.  Significance 
of  Defects  in  Welded  Structures,  T.  Kanayawa  and  A.  S.  Kobayashi,  ed..  University  of  Tokyo  Press,  1974,  p.  127-143. 

3.  SHAH,  R.  C.  Stress  Intensity  Factors  for  Through  and  Part-Through  Cracks  Originating  at  Fastener  Holes.  Mechanics  of  Crack 
Growth,  J.  R.  Rice  and  P.  C.  Paris,  ed.,  ASTM  STP  590,  1976,  p.  429-459. 

4.  CARTWRIGHT,  D.  J.  Stress  Intensity  Factors  and  Residual  Static  Strength  in  Certain  Structural  Elements.  Ph.D.  Thesis,  Depart- 
ment of  Mechanical  Engineering,  University  of  Southampton,  1971. 

5.  TIMOSHENKO,  S.,  and  GOODIER,  J.  N.  Theory  of  Elasticity.  McGraw-Hill  Book  Company,  Inc.,  2nd  Edition,  1951. 


Substitution  of  Equation  3 into  Equation  2 and  integrating  between  the  appropriate 
limits  gives: 


[ives;  _ \2  1 1 

_ go  r/R-d/2)(  177)  1*  (r7F/  j 

(r/R-l) 


(4) 


By  allowing  r ■ a + R to  be  substituted  in  Equation  4 and  equating*  which  in  turn 
is  equivalent  to  Equation  1 results  in: 


Kj/oq/ttsT 


(t^7r)  I (hot)2. 


(5) 


a/R 


In  Reference  1 the  authors  examine  the  biaxial ly-stressed  infinite  plate  with 
two  cracks  located  at  the  hole  boundary  as  shown  in  Figure  lb.  The  elastic  stresses 
in  the  absence  of  the  cracks  are  given  by:5 


o(r)  3 o0[l+(R/r)2]; 

the  mean  stress  over  the  crack  length  a is 
gp (a/R+2)  . 


(6) 


o_  = 
m 


la/KIT”’ 

and  hence,  proceeding  as  before,  the  normalized  stress  intensity  factor  is 


Kj/oq^  - 


(8) 


T> 


Figure  1 . Infinite  plate  contain- 
ing radial  cracks  at  the  boundary 
of  an  internal  circular  hole. 


°b 


•The  plane  stress  condition  is  cr  isidered  here. 


8 
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mM 


To  determine  the  accuracy  of  these  solutions  the  normalized  stress  intensity 
factor  resulting  from  Equations  5 and  8 are  compared  to  that  obtained  by  Bowie,6 
who  solved  the  same  cases  by  the  more  exact  method  of  conformal  mapping.  Such  a 
comparison,  showing  percent  differences  for  values  of  a/R  from  0 to  °°,  are  shown 
in  Table  1. 

It  is  seen  from  the  results  of  Table  1 that  for  both  cases  at  short  crack 
lengths,  a/R  = 0,  the  error  is  -12%;  and  as  a/R  increases  the  error  becomes  in- 
creasingly positive  up  to  a maximum  and  then  reduces  to  zero  when  a/R  = ■».  Compen- 
sation for  an  edge  crack  when  a/R  is  small  could  be  provided,  since  it  is  known' 
that  for  such  a configuration  normalized  stress  intensity  factor  approaches  1,12. 
This  will  however  increase  the  overestimate  at  larger  crack  ratios.  Nevertheless, 
the  method  provided  by  Williams  and  lsherwood1  appears  to  provide  estimates  of 
stress  intensity  factors  within  engineering  accuracy,  i.e.,  ±20%. 

In  Reference  1 the  authors  solve  several  other  problems  of  interest,  such  as 
the  determination  of  stress  intensity  factors  for  an  edge  crack  in  a finite  width 
plate  due  to  bending,  an  eccentrically-loaded  plate  with  unequal  cracks,  single- 
euge-notched,  center-notched,  double-edgc-notched  plates,  and  an  eccentric  center 
crack,  as  well  as  a radial  crack  at  the  bore  of  rotating  disk.  Further  details  of 
such  analyses  can  be  obtained  from  Reference  1. 

The  method  utilized  in  Reference  2 will  now  be  demonstrated.  Consider  Fig- 
ure 2a,  a uniaxial ly  stressed  infinite  plate  with  two  radial  cracks  at  the  boundary 
of  the  hole.  The  stress  intensity  factor  associated  with  this  figure  is  approxi- 
mately equivalent  to  that  shown  in  Figure  2b.  The  stress  distribution  Oyy(.x,0) 
in  an  uncracked  plate  with  a hole  of  radius  R is  determined  at  the  iocation  of  the 
cracks  for  the  applied  loading  in  question.  The  crack  of  length  2a,  as  shown  in 


Table  1.  COMPARISON  OF  APPROXIMATE  RESULTS  TO  MORE  EXACT  DATA 
FOR  AN  INFINITE  PLATE  CONTAINING  RADIAL  CRACKS 

Kj/oq/*# 


Two  Radial  Cracks 


a/R 

Uniaxial  Stress 

Biaxial  Stress 

Ref.  6 

Eq.  5 

Oiff . , X 

Eq.  11 

Oiff.,  % 

Ref.  6 

Eq.  8 

Oiff.,  %* 

0 

3.39 

3.00 

-12 

2.70 

-10 

2.26 

2.0 

-12 

0.10 

2.73 

2.70 

-1 

2.37 

-13 

1.98 

1 .91 

-4 

.20 

2.41 

2.47 

+2 

2.13 

-12 

1 .83 

1 .83 

0 

.30 

2.15 

2.29 

+7 

1.96 

-9 

1.70 

1.77 

+4 

.40 

1.96 

2.15 

+ 10 

1.82 

-7 

1.61 

1.71 

+6 

.50 

1.83 

2.04 

+ 11 

1.72 

-6 

1 .57 

1.67 

+6 

.60 

1.71 

1.94 

+13 

1.64 

-4 

1.52 

1.63 

+ 7 

.80 

1.58 

1.80 

+ 14 

1.51 

-4 

1 .43 

1.56 

+9 

1.00 

1.45 

1 .69 

+ 17 

1.43 

-1 

1.38 

1.50 

+9 

1.5 

1.29 

1.51 

+ 17 

1.28 

-1 

1.26 

1.40 

+ 11 

2.0 

1.21 

1.41 

+17 

1.23 

+2 

1 .20 

1 .33 

+11 

3.0 

1.14 

1.29 

+ 13 

1.16 

+2 

1.13 

1.25 

♦ 11 

5.0 

1.07 

1.18 

+ 10 

1.10 

+3 

1 .06 

1.17 

+ 10 

10.0 

1.03 

1.10 

+7 

1.05 

+2 

1.03 

1.09 

+6 

OD 

1.00 

1.00 

0 

1.00 

0 

1.00 

1.00 

0 

6.  BOWII-  , O.  L.  Analysis  of  an  Infinite  llate  Containing  Radial  l>acks  Originating  at  the  Boundary  of  an  Internal  Circular  Hole. 

I.  Math,  and  Phys..  v.  XXXV,  no.  II,  1956.  p.  60-71. 

7.  PARIS,  P.  C.,  and  Sill,  G.  C.  Stress  Analysis  of  l>acks.  Fracture  Toughness  Testing  and  Its  Applications.  ASTM  STP  381.  1970, 
p.  30-81. 


In  the  development  of  this  method  the  crack  is  assumed  to  be  absent  and  re- 
sidual surface  tractions  on  the  crack  surfaces  are  eliminated.  However,  in  deriv- 
ing Equation  5 the  mean  stress  arising  because  of  the  hole  is  applied  to  only  one 
radial  crack,  whereas  in  arriving  at  Equation  11  the  stress  due  to  the  hole  is 
integrated  across  both  crack  locations.  (Reference  3 allows  one  to  convert  from 
one  radial  crack  to  two  radial  cracks  at  the  hole  boundary  by  the  following 
relationship: 


one 

crack 


^ (l+a/2R)/ 1+a/R 


K 


• ) 

two 

cracks 


Thus,  the  method  of  Kobayashi2  and  Shah,3  although  more  difficult  to  apply,  is 
more  realistic  than  that  of  Williams  and  Isherwood.1  Nevertheless,  both  approaches 
will  provide  sufficient  accuracy  for  engineering  applications. 


The  approach  suggested  by  Cartwright1*  and  adapted  by  Baratta8  will  now  be 
outlined.  The  procedure  is  one  of  applying  the  appropriate  stress,  associated 
with  the  defect  discontinuity,  to  the  crack  tip  for  the  known  limiting  cases  of 
stress  intensity  and  utilizing  an  interpolation  scheme  between  limits  to  determine 
intermediate  values.  It  was  shown  by  the  preceding  illustrative  problems  how  the 
stress  intensity  could  be  approximated  by  applying  to  the  crack  faces  the  stresses 
that  would  normally  exist  on  the  plane  of  the  crack,  if  the  crack  was  absent. 
However,  in  the  analysis  to  follow,  rather  than  apply  the  stress  distribution  to 
the  crack  face,  only  the  stress  caused  by  the  defect  discontinuity  will  be  applied 
to  the  crack  tip. 


To  demonstrate  the  details  of  the  procedure  and  evaluate  the  validity  of  the 
approach,  it  will  be  applied  to  a problem  already  solved  by  more  accurate  means. 
Consider  a crack  in  a region  of  extreme  stress  concentration,  such  as  a single 
radial  crack  originating  from  a circular  hole  in  an  infinite  plate  subjected  to 
uniaxial  uniform  tension  at  infinity  as  shown  in  Figure  3a.  This  problem,  pre- 
viously used  as  an  example,  will  be  used  to  evaluate  the  following  approach. 


We  initially  consider  the  crack  length  a to  be  very  small  compared  to  the 
hole  radius  R.  Since  the  hole  radius  is  very  large  compared  to  the  crack  length 
we  can  consider  it  to  be  at  an  edge  of  a semi-infinite  plate,  as  shown  by  the 
equivalent  system  in  Figure  3b,  but  with  a stress  distribution  o0(r)  that  arises 
because  of  the  hole.  The  normalized  stress  intensity  ratio  for  an  edge  crack  in 
a semi-infinite  plate  stressed  in  tension  at  infinity  given  by  Paris  and  Sih7  is 
Kj/ao/im  = 1.12,  where  oq  is  the  remotely  applied  uniform  tensile  stress.  At  the 
opposite  extreme,  when  the  crack  length  is  very  large  compared  to  the  hole  radius, 
i.e.,  a/R  -*  <*>,  the  normalized  stress  intensity  ratio  Kj/ot/iTa  = 0.707,  from  Refer- 
ence 6.  It  appears  that  this  ratio  is  dependent  upon  some  function  of  a/R  which 
has  limits  of  1.12  when  a/R  = 0 and  0.707  when  a/R  ■>  °°.  Thus,  intermediate  values 
of  the  normalized  stress  intensity  ratio  can  be  approximated  by  choosing  a function 
which  is  well  behaved  between  the  above  stated  limits.  The  function  arbitrarily 
chosen  for  this  problem  as  well  as  those  to  follow  is  of  the  form: 


8.  BARATTA,  F.  I.  Stress  Intensity  Factor  Estimates  of  a Peripheral!  vOackrd  Spherical  Void  and  a Hemispherical  Surface  Pit. 
).  Am.  Ceram.  Soc.,  v.  61 . no.  1 1 , 1978,  p.  490-493;  also  Atmy  Materials  and  Mechanics  Research  ('enter,  AMMRC  TR  78-5, 
March  1978. 
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°o 


°b 

j Cracked  Circular  Hole 


ti.  Iquivalent  Cracked  Circular  Hole 


Figure  3.  Circular  hole  with  a 
single  radial  crack. 


f(a/R)  = Q - M(tan  'a/R)1”.  (12) 

Thus,  for  a circular  hole  with  one  radial  crack  the  normalized  stress  inten- 
sity ratio  is: 

Kj/oCr)*^  ■ f(a/R)  = Q-M(tan-1a/R)m  (13) 

where  Og(r)  is  the  stress  distribution  caused  by  the  circular  hole,  provided  by 
Timoshenko  and  Goodier5  which  is  Equation  3 in  this  paper.  Substitution  of  Equa- 
tion 3 into  Equation  13  and  transforming  the  coordinate  axis  r by  R + a,  gives 

Kj/oq^I  - [Q-M(tan-1a/R)m]  [ 1 + 0/2)  ( ) *+(3/2)  ( )“  ] . (14) 

Equation  14  represents  the  estimated  normalized  stress  intensity  ratio  Kj/o/iuT 
for  a radial  crack  of  length  a emanating  from  a circular  hole  of  radius  R.  The 
particular  constants  which  appear  in  this  equation  are  readily  determined  by  know- 
ing Q « 1.12  when  a/R  « 0;  M and  m are  found  by  allowing  f(a/R)  = 0.707  when 
a/R  * »>  and  f(a/R)  * 0.94  when  a/R  ■ 3.0  (taken  from  Reference  6) . With  f(a/R)  = 
0.94,  the  error  span  on  the  function  would  be  minimized  throughout  the  total  range 
of  a/R.  Thus,  Q ■ 1.12;  M was  found  to  be  0.119  and  m equals  2.748. 

Equation  14  was  then  utilized  to  obtain  K/o/tuT  as  a function  of  a/R  from  0 
to  ®.  It  is  seen  that  these  results,  shown  in  Table  2,  are  within  ±0%  when  com- 
pared to  the  more  accurate  data  of  Reference  6.  Although  the  results  of  this  pro- 
cedure are  approximate,  it  is  believed  that  the  errors  involved  when  applied  to 
other  analogous  problems  will  be  of  similar  magnitude  as  those  given  in  the  table. 
Thus,  a similar  approach  is  used  to  solve  these  problems  of  interest. 
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Application  of  the  above  technique 
will  be  demonstrated  by  providing  esti- 
mates of  stress  intensity  factors  for  a 
spherical  void  with  a circumferential 
crack  at  its  equator  in  an  infinite  body 
under  uniaxial  tensile  stress,  as  shown 
in  Figure  4.  Also  included  is  a circum- 
ferentially-cracked hemispherical  pit  at 
a free  surface  of  a semi-infinite  body 
also  stressed  by  uniaxial  tension,  shown 
in  Figure  5.  Since  these  types  of  de- 
fects occur  in  both  metallic  and 
ceramic8  systems  such  solutions  are 
of  interest. 


Table  2.  COMPARISON  OF  APPROXIMATE  TO  EXACT  STRESS 
INTENSITY  RESULTS  FOR  A CIRCULAR  HOLE 
WITH  RADIAL  CRACK 


Kj/o/itT 


Goodier9  provided  the  stress  dis- 
tribution for  a spherical  defect  in  an 
infinite  medium  stressed  in  tension. 
Because  of  the  defect,  three  stresses 
are  developed,  which  are  indicated 


a/R 

From 
Eq.  14 

Ref.  6 

Diff.,  % 

0 

3.36 

3.39 

-1 

0.1 

2.73 

2.73 

0 

.2 

2.32 

2.30 

+1 

.3 

2.04 

2.04 

0 

.4 

1.83 

1.86 

-2 

.5 

1.69 

1.73 

-2 

.6 

1.56 

1.64 

-5 

.8 

1 .40 

1.47 

-5 

1.0 

1.29 

1.37 

-6 

1.5 

1.13 

1.18 

-4 

2.0 

1.03 

1.06 

-3 

3.0 

0.94 

0.94 

0 

5.0 

0.86 

0.81 

+6 

10.0 

0.78 

0.75 

+4 

o° 

0.708 

0.707 

0 

Figure  S.  A hemispherical  surface  pit  with  a 
circumferential  crack  at  its  semi-equator. 


Figure  4.  A spherical  void  with  a 
circumferential  crack  at  its  equator. 


9.  GOODIER,  J.  N.  Concentration  of  Stress  Around  Spherical  and  Cylindrical  Inclusions  and  Flaws.  Trans.  ASME,  v.  55,  1933, 

p.  39-44. 
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schematically  in  Figure  6.  They  are:  a radial  stress  or  acting  normal  to  the 
interface,  a stress  o0  acting  in  a tangential  direction  to  a meridian  (the  north- 
south  pole  axis  of  the  sphere  is  considered  aligned  with  the  axis  of  the  applied 
stress),  and  another  tangential  stress  normal  to  both  or  and  oe.  In  the  follow- 
ing analysis,  the  maximum  tensile  stress  is  assumed  to  cause  crack  initiation. 

It  is  indicated  in  Reference  9 that  the  maximum  tensile  stress  for  a spher- 
ical void  oacurs  when  6*0,  resulting  in  a "hoop  stress"  o0  which  girdles  the 
equator.  It  is  expected  that  this  tensile  stress,  if  large  enough,  can  initiate 
an  axisymmetric  crack  and  extend  it  in  the  radial  direction  as  envisioned  in  Fig- 
ure 4.  The  stress  distribution  when  6 = 0 is  given  in  the  following  equation: 

vr>  ■ °o  j [2Ttw]  Msnn’H-  (is) 

where  v is  Poisson's  ratio  of  the  material  and  r is  the  radial  distance  from  the 
center  of  the  spherical  void  (r  > R). 


To  obtain  an  approximation  for  the  stress  intensity  factor  associated  with 
Figure  4,  we  again  initially  assume  that  the  crack  length  a is  small  compared  to 
the  sphere  radius  R.  Also  as  before,  we  further  assume  that  when  the  crack  is 
small,  the  influence  of  the  spherical  void  on  the  crack  is  negligible  and  thus 

can  be  considered  as  an  edge  crack  as  shown  in  Figure  3b. 
Now  the  stress  distribution  o0(r)  is  provided  by  Equa- 
<r0  tion  15,  which  is  due  to  the  spherical  void.  Also,  we 

assume  that  Equation  12  is  applicable  as  well  as  the  con- 
stants (except  M)  previously  given  as  Q = 1.12  and  m * 

<rr  2.748.  M is  determined  by  realizing  that  in  the  extreme 


limit  when  a/R  “>  the  defect  geometry  approaches  a disk 
crack  and  thus  Kj/o/rTa  -*•  2/tr,  according  to  Sneddon.10 
Therefore,  making  use  of  this  end  limit  in  Equation  12, 
which  is  substituted  along  with  Equation  15  into  Equa- 
tion 13,  results  in  the  following  normalized  stress  inten- 
sity expression: 


’ where  Q = 1.12,  m = 2.748,  and  M is  determined  to  be  0.101. 

o 

Figure  6.  Coordinate  system  to 
describe  the  stress  state  around 
a spherical  void. 


10.  SNEDDON,  I.  N.  The  Distribution  of  Stress  In  the  \etghborht>od  of  o Crock  in  on  Has  tic  Solid.  Proc.  Roy.  Soc.  London,  v.  A- 1 87, 
1946. 
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The  resulting  calculations  using  Equation  16,  which  represents  the  normalized 
stress  intensity  ratio  for  a spherical  void  cracked  at  its  equator  and  stressed  by 
uniform  uniaxial  tension  at  infinity,  are  shown  in  Table  3 as  a function  of  a/R 
from  0 to  «>  and  for  v = 0.25  and  0.30. 

The  final  case  considered  is  a hemispherical  pit  at  a free  surface  of  a semi- 
infinite body  with  a crack  initiated  at  the  semi-equator  of  the  pit  and  extended 
radially  by  a uniaxial  uniform  tensile  stress  as  shown  in  Figure  5.  As  outlined 
previously,  we  initially  assume  that  a is  very  small  compared  to  R.  As  before,  we 
further  assume  that  when  the  crack  is  small,  the  influence  of  the  hemispherical 
pit  on  the  crack  is  negligible  and  thus  can  be  considered  an  edge  crack  as  shown 
in  Figure  3b,  with  the  stress  distribution  ae(r)  caused  by  the  surface  pit  discon- 
tinuity. Again  as  before,  we  utilize  Equation  13,  with  Q = 1.12  and  m = 2.748, 
but  determine  M from  the  end  limit  when  a/R  -+•  °°.  Smith  and  Alavi11  give  the  varia- 
tion of  Kj/o/im  around  the  semicircular  crack  boundary  and  indicate  that  the  maxi- 
mum value  occurs  at  the  junction  of  the  crack  with  the  free  surface  and  is  1.21  x 
2/tt . This  results  in  M = 0.140. 

The  stress  distribution  ae(r)  is  provided  by  Eubanks,12  who  determined  the 
variation  of  the  circumferential  stress  along  the  axis  of  symmetry  as  a function 
of  r at  the  base  of  the  pit.  Equation  12,  including  the  constants  previously  given, 
i.e.,  Q = 1.12,  m = 2.748,  and  M = 0.140,  and  the  normalized  stress  distribution 
o0(r)/oo  given  in  tabular  form,  will  provide  the  normalized  stress  intensity  ratio 
for  a peripherally-cracked  hemispherical  surface  pit  with  v = 0.25.  Such  results, 
including  ae(r)/oo,  are  given  in  Table  4 as  a function  of  a/R  from  0 to  «. 


The  application  of  Equation  14  and  Tables  3 and  4 has  certain  restrictions 
resulting  from  the  implied  assumptions  used  in  References  9 and  12  when  determin- 
ing the  appropriate  stress  distribution.  These  restrictions  are: 

The  material  is  isotropic  and  homogenous. 

The  defect,  which  includes  the  crack,  is  very  large  compared  to  the 
grain  size  of  the  crystalline  aggregate. 


Table  3.  STRESS  INTENSITY  RATIO  FOR 
A PERIPHERALLY-CRACKED  SPHERICAL  VOID 


Kj/jo/jrT 


a/fK 

v 0.25 

0.30 

0 

2.26 

2.30 

0.15 

1.73 

1.75 

0.35 

1.41 

1.41 

0.55 

1 .26 

1.26 

1.00 

1.10 

1.08 

2.00 

0.94 

0.93 

3.00 

.86 

.86 

5.00 

.79 

.79 

10.00 

.72 

.72 

00 

.64 

.64 

Table  4.  NORMALIZED  STRESS  INTENSITY  RATIO 
FOR  A PERIPHERALLY -CRACKED  HEMISPHERICAL 
SURFACE  PIT,  WITH  v = 0.25 


a/R 

Kj /ofrd 

0 

1.23 

2.50 

0.15 

1.63 

1.82 

0.35 

1.29 

1.44 

0.55 

1.15 

1.27 

1 .00 

1.04 

1.11 

2.00 

1 .00 

0.99 

3.00 

1.00 

.93 

5.00 

1.00 

.88 

10.00 

1.00 

.83 

00 

1.00 

.77 

11.  SMITH,  F.  W„  and  ALAVI,  M.  1.  Stress  Intensity  Factors  for  a Part-Through  Circular  Surface  Flaw.  Pioc.  1st  Inti.  Conf.  on 
Pressure  Vessel  Tech.,  Delvt,  Holland.  1969. 

12.  EUBANKS,  R.  A.  Stress  Concentration  Due  to  Hemispherical  Pit  at  a Free  Surface.  J.  Appl.  Mcch.,  v.  21,  19S4,  p.  57-62. 
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The  tip  of  the  crack,  which  surrounds  the  spherical  void,  can  be  no 
closer  than  2R  + 2a  to  a body  boundary. 

The  tip  of  a crack,  which  surrounds  a hemispherical  pit  at  a free 
surface,  can  be  no  closer  than  R + a away  from  any  other  body  boundary. 

Since  there  are  no  other  available  solutions  to  effect  a comparison  of  the 
three-dimensional  cases  given  above  their  accuracy  can  not  be  definitively  ascer- 
tained. However,  the  author0  of  this  technique  claims  that  engineering  accuracy 
is  preserved. 

SUMMARY 

It  is  seen  that  the  method  used  by  Williams  and  Isherwood1  as  applied  to  two 
radial  cracks  emanating  from  a circular  hole  stressed  by  uniaxial  tension  resulted 
in  a range  of  error  of  -12%  to  +17%.  For  this  particular  class  of  problems  the 
approach  does  not  allow  differentiation  between  the  one-crack  case  and  two-crack 
case.  This  method  should  be  restricted  in  application  to  those  problems  where  only 
one  crack  is  present. 

The  more  complex  method  of  Kobayashi2  and  Shah3  will  allow  such  flexibility 
and  is  more  accurate,  i.e.,  the  error  ranged  from  -13%  to  +3%  for  the  same  prob- 
lem. However,  for  the  two-crack  biaxial  stress  case  the  range  of  error,  when 
compared  to  Bowie's  results6  was  -12%  to  +11%.  Nevertheless,  both  methods  will 
provide  results  that  can  be  adapted  for  engineering  applications. 

The  approach  suggested  by  Cartwright1*  and  adapted  by  Baratta0  is  even  easier 
to  apply  but  one  must  know  the  end  limit  to  obtain  an  interpolation  relationship. 
The  resulting  accuracy  is  not  definitively  known  for  the  three-dimensional  cases 
solved  but  it  is  believed  that  it  should  be  well  within  acceptable  standards 
for  engineering  applications. 
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